Solar Hydrogen Burning and Neutrinos 



W.C. Haxton\ P.D. Parker^, C.E. Rolfs^ 



^ INT and Physics Department, Univ. Washington, Seattle, WA, 98195-1550 USA 
^ Physics Department, Yale University, New Haven, CT, 06520-8124 USA 
l/^ , ^ Experimentalphysik III, Ruhr- Universitdt Bochum, D-44'^^0 Bochum, Germany 

o 
o 

(N 

I Abstract 

We summarize the current status of laboratory measurements of nuclear cross sec- 
tions of the pp chain and CN cycle. We discuss the connections between such 
measurements, predictions of solar neutrino fluxes, and the conclusion that solar 
I neutrinos oscillate before reaching earth. 
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^ . 1 Introduction 

The probable role of thermonuclear reactions in generating energy within the 
^ . solar interior was recognized as early as the 1920s by astrophysicists such as 

Eddington [1]. This deduction came about by default, rather than by any direct 
laboratory observation of the reactions. Once the earth's age was established 
to be at least a few billion years old (see Burchfield [2] for an interesting history 
of the development of this idea), it became apparent that gravitational and 
chemical energy sources could not sustain solar evolution, given the known 
mass, radius, and luminosity of the sun. In contrast, newly discovered nuclear 
reactions and decays could readily supply the required energy over the requisite 
time. A key step was taken by Gamow [3], who recognized that tunnelling 
under the Coulomb barrier could enable nuclear reactions to occur at the low 
energies typical of the solar interior, where kT ~ a few keV {e.g., Atkinson & 
Houtermans [4]). This idea and subsequent investigations led to descriptions of 
possible reaction chains for stellar energy generation. By the end of the 1930s 
Weizsacker [5,6], Bethe & Critchfield [7], and Bethe [8] had described the 
essential features of the carbon-nitrogen cycle and the proton-proton chain, 
the processes by which the sun and similar stars convert hydrogen to helium. 
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The neutrino has a similar vintage, and was connected early on to stellar 
reactions. In December 1930 Wolfgang Pauli proposed that the emission of 
an unobserved spin-1/2 neutral particle - the electron neutrino (ve) ~ might 
explain the apparent lack of energy conservation in nuclear beta decay 

{A, Z) ^ {A,Z -l) + e+ + ue 

Enrico Fermi was present at a number of Pauli 's presentations and discussed 
the neutrino with him on these occasions. In 1934, following closely Chad- 
wick's discovery of the neutron, Fermi proposed a theory of beta decay based 
on Dirac's description of electromagnetic interactions, but with weak currents 
interacting at a point, rather than at long distance through the electromag- 
netic field. He described beta decay as a proton decaying to a neutron, a 
phenomenon energetically possible because of nuclear binding energies, with 
the emission of a positron and neutrino. Apart from the absence of parity 
violation, which would not be discovered until 1957, Fermi's description is 
a correct low-energy approximation to our current standard model of weak 
interactions. 

Thus, as nuclear astrophysicists unravelled the stellar processes for hydrogen 
burning, stars were recognized to be copious sources of neutrinos 

4p ^ ^He + 2e+ + 

This observation led Ray Davis to construct the chlorine experiment to mea- 
sure the flux of solar neutrinos. The resulting solar neutrino puzzle was ul- 
timately traced to new physics (neutrino oscillations) requiring the mixing 
of massive neutrinos, a phenomenon beyond the standard electroweak model. 
This discovery depended critically on the reliability of solar model neutrino 
flux predictions, and thus on the quality of our understanding of the nuclear 
reaction rates governing the pp and CN cycles. 

The goal of this paper is to summarize our understanding of solar hydro- 
gen burning, with special emphasis on the connections between laboratory 
astrophysics — measurement of nuclear cross sections at or near energies char- 
acteristic of the solar core — and the predictions of solar neutrino fluxes. 



2 Hydrogen Burning 

Observations of stars reveal a wide variety of stellar conditions, with lumi- 
nosities relative to solar spanning a range L ~ 10~^ to 10^ Lq and surface 
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temperatures Tg ■--'2000-50000 K. The simplest relation one could propose 
between luminosity L and is that for a blackbody 

L = AirR'aT^ L/L^ = (R/Re)' {Ts/T^f (1) 

which suggests that stars of a similar structure might lie along a one-parameter 
path (corresponding to R/ Rq above) in the luminosity (or magnitude) vs. tem- 
perature (or color) plane. In fact, there is a dominant path in the Hertzsprung- 
Russell color-magnitude diagram along which roughly 80% of the stars reside. 
This is the main sequence, those stars supporting themselves by hydrogen 
burning through the pp chain or CN cycles. 

As one such star, the sun is an important test of our theory of main sequence 
stellar evolution: its properties - age, mass, surface composition, luminosity, 
and helioseismology - are by far the most accurately known among the stars. 
The standard solar model (SSM) traces the evolution of the sun over the past 
4.6 billion years of main sequence burning, thereby predicting the present- 
day temperature and composition profiles, the relative strengths of competing 
nuclear reaction chains, and the neutrino fluxes resulting from those chains. 
Standard solar models share four basic assumptions: 

• The sun evolves in hydrostatic equilibrium, maintaining a local balance 
between the gravitational force and the pressure gradient. To describe this 
condition in detail, one must specify the equation of state as a function of 
temperature, density, and composition. 

• Energy is transported by radiation and convection. While the solar enve- 
lope is convective, radiative transport dominates in the core region where 
thermonuclear reactions take place. The opacity depends sensitively on the 
solar composition, particularly the abundances of heavier elements. 

• Thermonuclear reaction chains generate solar energy. The standard model 
predicts that over 98% of this energy is produced from the pp chain con- 
version of four protons into ^He (see Fig. 1), with proton burning through 
the CN cycle contributing the remaining 2% (Fig. 2). The Sun is a large 
but slow reactor: the core temperature, T^. ~1.5 x 10^ K, results in typical 
center-of-mass energies for reacting particles of ~ 10 keV, much less than 
the Coulomb barriers inhibiting charged particle nuclear reactions. Thus 
reaction cross sections are small: in most cases, laboratory measurements 
are only possible at higher energies, so that cross section data must be 
extrapolated to the solar energies of interest. 

• The model is constrained to produce today's solar radius, mass, and lu- 
minosity. An important assumption of the standard model is that the Sun 
was highly convective, and therefore uniform in composition, when it first 
entered the main sequence. It is furthermore assumed that the surface abun- 
dances of metals (nuclei with A > 5) were undisturbed by the subsequent 
evolution, and thus provide a record of the initial solar metalicity. The re- 
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maining parameter is the ^He/H ratio, which is adjusted until the model 
reproduces the present solar luminosity after 4.6 billion years of evolution. 
The resulting ^He/H mass fraction ratio is typically 0.27 ± 0.01, which can 
be compared to the big-bang value of 0.23 ± 0.01. Note that the Sun was 
formed from previously processed material. 

The model that emerges is an evolving Sun. As the core's chemical composition 
changes, the opacity and core temperature rise, producing a 44% luminosity 
increase since its arrival on the main sequence. The temperature rise governs 
the competition between the three cycles of the pp chain: the ppl cycle domi- 
nates below about 1.6 x 10^ K; the ppll cycle between (1.7-2.3) x 10^ K; and 
the ppIII above 2.4 x 10^ K. The central core temperature of today is about 
1.55 X 10^ K. 

This competition between the cycles determines the pattern of neutrino fluxes. 
One consequence of the thermal evolution of our sun is the relatively recent 
appearance of a significant neutrino flux, the species that dominated the 
Davis experiment. As this flux is produced in the ppIII cycle, it is a very 
sensitive function of the solar core. The SSM predicts an exponential increase 
in this flux with a doubling period of about 0.9 billion years. 

A final aspect of SSM evolution is the formation of composition gradients on 
nuclear burning timescales. Clearly there is a gradual enrichment of the solar 
core in ^He, the ashes of the pp-chain. Another element, ^He, is produced and 
then consumed in the pp chain, eventually reaching some equilibrium abun- 
dance. The time-scale for equihbrium to be estabhshed as well as the eventual 
equilibrium abundance arc both sharply decreasing functions of temperature, 
and thus increasing functions of the distance from the center of the core. Thus 
a steep ^He density gradient is established over time. 

The SSM has had some notable successes. From helioseismology the sound 
speed profile c(r) has been very accurately determined for the outer 90% of 
the Sun, and is in good agreement with the SSM. (However, this conclusion 
depends in part on inputs to the SSM: recent calculations employing one set 
of revised opacities do show some significant discrepancies, as we will discuss 
later.) Helioseismology probes important predictions of the SSM, such as the 
depth of the convective zone. The SSM is not a complete model in that it 
does not attempt to describe from first principles aspects such as convection 
that extend beyond ID descriptions. For this reason certain features of solar 
structure, such as the depletion of surface Li by two orders of magnitude, are 
not reproduced. Li destruction is often attributed to convective processes that 
operated at some epoch in our sun's history, dredging Li to a depth where 
burning takes place. 
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3 Reaction Rate Formalism 



3.1 Nuclear rates and S- factors 



Perhaps the key step that allowed quantitative modelling of the sun and 
other main-sequence stars was the description of the relevant nuclear reac- 
tion chains. This required both the development of a theory for sub-barrier 
fusion reactions and careful laboratory measurements to constraint the rates 
of these reactions at stellar temperatures. 

At the temperatures and densities in the solar interior {e.g., Tc ~ 15 x 10^ 

K and pc ~ 153 g/cnv' [9]), interacting nuclei reach a Maxwellian equilibrium 
distribution in a time that is infinitesimal compared to nuclear reaction time 
scales. Therefore, the reaction rate between two nuclei can written as {e.g., 
Burbidge et al, [10]; Clayton [11]; Rolfs & Rodney [12]): 

nin2 , > . s 

ri2 = T— Wl2- (2) 
1 + 0i2 

Here the Kronecker delta prevents double counting in the case of identical 
particles, rii and n2 are the number densities of nuclei of type 1 and type 2 
(with atomic number Zi and Z2, and atomic mass Ai and A^), and {o'v)^2 
denotes the product of the reaction cross-section and the relative velocity v 
of the interacting nuclei, averaged over the collisions in the stellar gas, 

00 

{av)u = j[(7{v) ■ v\^{v)dv. (3) 



Under solar conditions nuclear velocities are very well approximated by a 
Maxwell-Boltzmann distribution. It follows that the relative velocity distri- 
bution is also a Maxwell-Boltzmann, governed by the reduced mass of the 
colliding nuclei. 



Therefore, 

where E is the relative kinetic energy in the center-of-mass system. In order 
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to evaluate {av) the energy dependence of the reaction cross-section must be 
determined. (In Fowler et al. [13] and more recently in Angulo et al. [14]. the 
appropriate expressions for {av) are tabulated for a variety of specific reac- 
tions, including both non-resonant and resonant neutron-induced and charged- 
particle-induced reactions.) 

Almost all of the nuclear reactions relevant to solar energy generation are 
nonresonant and charged-particle induced. For such reactions, as the energy 
dependence at low energies is dominated by the Coulomb barrier, cross sections 
are typically expressed in terms of the S-f actor, 

^{E) exp (-27r77) , (6) 



where the 1/E term comes from the geometrical factor, exp(— 27rr7) is the 
usual Gamow barrier penetration term with 77 = Z1Z2 a/v, v is the relative 
velocity, and a ~ 1/137 is the fine structure constant. (We take h — c — 
1.) That is, the sharp energy dependence associated with s-wave interactions 
of point nuclei has been removed from ct{E), leaving a quantity S{E) that 
contains the nontrivial nuclear physics, yet is simpler to model over an energy 
range spanning both solar reactions and terrestrial laboratory measurements. 
For nonresonant reactions S{E) will be a slowly varying function of E. 

One must take into account differences in the atomic environments to cor- 
rectly relate laboratory measurements of (^{E) and S{E) to the corresponding 
quantities in the solar interior. As light nuclei in the solar core are almost 
completely ionized, the solar electron screening correction /q can be treated 
in a weak-screening approximation [15]. The correction depends on the ratio 
of the Coulomb potential at the Debye radius Rd to the temperature, 

fo = exp (1^) = exp (0.188 Z, ( pT Te"^/^) , (7) 



where C,Rb — (if^j ) P is the number density of nucleons, po is a dimen- 

sionless density measured in g/cm^, C = {^i ^ + f^) ] -'^i is the 

mass fraction of nuclei of type i, and Tg is the dimensionless temperature in 
units of 10^ K. Note that the screening factor enhances solar cross sections. 



Substituting (4) into (1) yields 

1/2 00 



{av) 
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fo J S{E) exp {-27:7]) exp{-E/kT) dE. (8) 
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The energy dependence of this integrand is shown schematically in Fig. 3. In 
order to evaluate this integral we expand S{E) in a Taylor series, 

S{E) = S{0) + ES'{0) + -^^^"(O) + (9) 

Substituting this into (6) yields (Bahcall [16]), 

^"""^'^ = []Ikf) ^''^^^ ^'^^^ """^ i-^Eo/kT) (10) 
= 1.301 X 10-^^cmVs {^f'fo j^^^ T-^/^xp (-3^oAT) 
where 

Eo/kT ^ {nZ,Z2a/V2f'{picykTy/\ 
AE,/kT = A[Eo/3kTf\ A = 

and 

corresponds to the maximum of the integrand, the Gamow peak, and is 
thus the most probable energy of reacting nuclei (see Fig. 3). A£^o corresponds 
to the full width of the integrand at 1/e of its maximum value. (For example, 
for a typical solar core temperature of Tg = 16, Eq = 23.4 keV and AEo= 13.1 
keV for the ^He(a,7)''Be reaction in the proton-proton chain.) 

Numerous laboratory measurements have been carried out to determine S{0), 
S'{0), and S"'(0) for the various nuclear reactions contributing to the proton- 
proton chain and the CN cycle. At the present time it is not technically feasible 
to measure the cross sections for most of these reactions in the region of Eq, 
as the cross sections are severely suppressed by the Coulomb barrier. (Clearly, 
as the basic timescale is the solar age, the sun is a very slow reactor operating 
at a temperature where many reactions are highly suppressed.) Therefore Seff 
often must be determined from extrapolations from laboratory measurements 
at higher energies, typically E ^ 100 to 200 keV. 

So far, the one significant reaction for which it has been possible to measure 
the cross section down to Eq is the ^He(^He,2p)'^He reaction - see Section 4.3, 
below. At the sun's center, Eq ~ 22 keV for this reaction. The measured cross 
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section is a{Eo) ~1.5 pb. This can be compared with the cross sections for 
the (p,7) and {a,j) capture reactions that contribute in the pp chain or CN 
cycle, e.g., 

3He(Q;,7)^Be ^ Eo 23 keV ^ a(Eo) ~ 3 x 10"^ pb 

^Be(p,7)SB ^ - 18.4 keV ^ a{Eo) ~ 1.5 x lO^^ 

^^N(p,7)i50 ^Eor^ 27.2 keV ^ a(Eo) ~ 2.2 x 10"^ pb 

Direct measurements of these much smaller capture cross sections at E^ are 
not anticipated without orders of magnitude improvements in sensitivity. 

For completeness we can also evaluate Eq. (5) for the case of a reaction dom- 
inated by a narrow, isolated resonance, 

oo oo 

j Ea{E) exp {-E/kT) dE = E^es exp {-E^es/kT) j a {E) dE (11) 



rir2 



Assuming a Breit-Wigner line shape a (E) = (Jbw {E) = ttA^cj -2 — ^ 

with the proper statistical factor uj = {2j + l)/((2si + 1) (2s2 + 1)) yields 



00 

j gbw (E) dE 



UJ — — — = UJ 7. 

2 r 2 ' 



Expressing T in units of 10^ K and generalizing for a series of isolated reso- 
nances yields 

{av),, = 2.53 X 10-^=^cmVs {AT,)-"^ ^ Mi x exp {-E^/kT) . (12) 



For intermediate cases, between the limits of narrow, isolated resonances and 
non-resonant reactions, one must directly integrate Eq. (5). 



3.2 Tests of electron shielding 



In the above treatments, it is assumed that the Coulomb potential of the 
target nucleus and projectile is that resulting from bare nuclei. However, for 
nuclear reactions studied in the laboratory, the target nuclei and the projectiles 
are usually in the form of neutral atoms or molecules and ions, respectively. 
The electron clouds surrounding the interacting nuclides act as a screening 
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potential: the projectile effectively sees a reduced Coulomb barrier, both in 
height and radial extension. This, in turn, leads to a higher cross section for 
the screened nuclei, (Js{E), than would be the case for bare nuclei, ah{E). 
There is an enhancement factor (Assenbaum et al. [17]) 

fiah{E) = a,{E)/ab{E) ~ e-K^{T:r]U J E) > 1 for Ue « E, (13) 

where Ue is an electron-screening potential energy. This energy can be calcu- 
lated, for example, from the difference in atomic binding energies between the 
compound atom and the projectile plus target atoms of the entrance channel. 
Note that for a stellar plasma, the value of the bare cross section ab{E) must 
be known because the screening in the plasma will be quite different from that 
in the laboratory nuclear-reaction studies, i.e., (Tpiasma{E) = fpiasma{E) (Xb^E), 
where the plasma enhancement factor fpiasma{E) must be explicitly included 
for each situation. The later factor involves the Debye radius of the electrons 
around the ions in the plasma, as discussed earlier. Similarly, a good under- 
standing of electron-screening effects in the laboratory is needed to extract a 
reliable (Tb{E) from low-energy data taken with terrestrial atomic targets. A 
great deal can be done experimentally to test our understanding of electron 
screening. Thus, while nuclear reactions are not measured under the condi- 
tions found in stellar interiors, we can test of understanding of the screening 
corrections needed to determine ab{E) from laboratory data. 

Experimental studies of reactions involving light nuclides (Engstler et al. [18]; 
Strieder et al. [19] and references therein) have shown the expected exponential 
screening enhancement of the cross section at low energies. For example, the 
cross section of the ■^Hc(d,p)^He reaction was studied over a wide range of 
energies (Aliotta et al. [20]) where new energy loss information at the relevant 
low energies was available for the analysis of the data: the results led to Ue = 
219±15 eV, significantly larger than the adiabatic limit from atomic physics, 
Uad = 119 eV. A weak point in the analysis is the assumption of the energy 
dependence of the bare cross section ab{E). A direct measurement of (Tb{E) 
using bare nuclides {e.g., a crossed ion-beam set-up) appears difficult if not 
impossible due to luminosity problems. However, a new experimental method 
for an indirect (Jb{E) determination has been developed, the Trojan Horse 
Method (Strieder et al. [19] and references therein). Although it is not yet 
clear whether all relevant components of this method have been included in 
current analyses, the available data demonstrate how well the indirect method 
and the direct measurements complement one another. 

There exist various surrogate environments that have allowed experimentalists 

to test our understanding of plasma screening effects. Recently, screening in 
d(d,p)t has been studied for deuterated metals, insulators, and semiconduc- 
tors, i.e. 58 samples in total (Raiola et al. [21], Raiola et al. [22] and references 
therein). As compared to measurements performed with a gaseous D2 target 
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(Ue = 25 eV), a large effect has been observed in all metals (of order Ue — 300 
eV) , while a small (gaseous) effect is found for the insulators and semiconduc- 
tors. For the metals, the hydrogen-solubilities are small (a few percent) leaving 
the metaUic character of the samples essentially unchanged. An explanation of 
the large effects in metals was suggested by the classical plasma screening of 
Debye applied to the quasi-free metallic electrons. The electron Debye radius 
around the deuterons in the lattice is given by 

Rd = {sokT/aneSfPaf'' = 69(T/ne//Pa)'/' N (14) 

with the temperature of the free electrons T (in units of K) , ng/ / the number 
of valence electrons per metallic atom, and the atomic density Pa (in units 
of m~^). With the Coulomb energy between two deuterons Rd set equal 
to ?7e, one obtains Ue = (47r£o)~'^ a/Ro- A comparison of the calculated and 
observed Ue values determines ^^e//, which is for most metals of the order 
of one. The acceleration mechanism of the incident positive ions leading to 
the high observed Ue values is thus the Debye electron cloud at the rather 
small radius Ro, about one tenth of the Bohr radius. The rie// values have 
been compared with those deduced from the known Hall coefficient: within 
2 standard deviations the two quantities agree for most metals. A critical 
test of the Debye model is the predicted temperature dependence U^iT) oc 
T"^/^, which was also confirmed experimentally. Thus, metals appear to form 
a "plasma for the poor man," allowing terrestrial testing of screening effects 
very similar to those found in stellar plasmas. 

In the next sections we discuss the present status of the experimental mea- 
surements and extrapolations for each of the various nuclear reactions in the 
proton-proton chain and the CN cycle. 



4 The Proton— Proton Chain and The Carbon— Nitrogen Cycle 

4-1 The pp and pep reactions 

Unlike the typical stellar nuclear reactions for which rates are determined by 
extrapolating laboratory measurements to the needed Eq, the weak rates for 
the pp beta-decay and pep electron capture reactions are too small to be 
measured directly in the laboratory. (More correctly, the Sudbury Neutrino 
Observatory experimentalists measured the charged-current breakup reaction 
on the deuteron, the inverse of the pp reaction, which has been used to place 
some direct constraints on the size of the pp rate.) Instead the SSM values 
for these rates are based on model calculations that take into account the 
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measured values of GpcosOc and of the axial-vector coupling constant 
(determined from the neutron half-life and from superallowed 0^ beta- 
decay in nuclei) , the nuclear structure governing the formation of deuterium, 
exchange current contributions to the weak amplitude, and weak radiative 
corrections to that amplitude. 

One approach follows in spirit of early estimates of the pp rate {e.g., Salpeter 
[23] and Bahcall and May [24]), and is based on a potential-model descrip- 
tion of the deuteron. Modern treatments employ reahstic NN potentials that 
are derived from phase-shift analyses of scattering data. Perhaps the primary 
nuclear structure uncertainty comes from two-body exchange currents, which 
are dominated by short-ranged exchanges involving poorly known couplings. 
However these are tightly constrained by the known strength of tritium beta 
decay. As the nonrelativistic two- and three-nucleon problems can be solved 
exactly (and assuming there are no important three-body exchange currents), 
a very precise rate estimate can be obtained. 

More recently an approach has been taken based on pionless effective field 
theory that illustrates nicely the functional dependence of the pp rate on 
model-dependent short-range physics. In these studies [25,26,27] the S-factor 
is expressed as 

5„(0) = 3.93xlO-MeVb (^) (^)' (Sl^W (15) 
^ \6m J V 1-2670 J V 3073s ) 



where A is proportional to the hadronic axial matrix element connecting the 
pp and deuteron states, including the consequence of short-range two-nucleon 
physics. Effective field theory determines the functional form of A, 

A = 2.58 + 0.011 f^^) - 0.0003 f^^) (16) 
Ifm / Vlfm 



where Li^a and Ki a parameterize the two-nucleon physics at next-to-leading 
order (NLO) and next-to-next-to-next-to-leading order (NNNLO), respec- 
tively. As the correction due to K^^a is expected to be much less than 1%, 
it follows that the model dependence is embodied in Li ,4, a single parame- 
ter that determines both S(0) and its energy dependence. From the work of 
Schiavilla et al. [28], this is fixed as 4.2 ± 2.4 fm^ from tritium beta decay. It 
follows that 

Spp (0)=(3.92± 0.08) X 10-25 MeVb 

The central value is 1% smaller than that recommended by a 1997 INT working 
group [29]. The uncertainty encompasses that earlier result. 
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A small contribution to pp fusion comes from electron capture on corre- 
lated protons in the solar plasma. To the accuracy required, the rate for 
p + e^ + p ^ d + i^e is proportional to that for pp beta-decay, Rpp. Bahcall 
and May [24] found 

Rpep = 5.51 X 10-' p (1 + X^)T6"'/'(1 + 0.02T6) Rpp (17) 



where Xh is the mass fraction of hydrogen. In the solar interior Rpep ~ 2.4 x 

10~^ Rrm- 



4-2 The d(p, ^)^He reaction 



Early measurements (Griffiths et al. [30] and Schmid eA, al. [31]), utilizing D2O 
ice targets sufficiently thick to stop the incident protons, were used to extract 
d(p,7) cross sections down to center-of-mass energies as low as 10 keV. More 
recently, using a differentially pumped gas target in an underground laboratory 
(LUNA — See Section 4-3)-, Casella et al. [32] have measured the D(p,7)^He 
reaction down to Ecm = 2.5 keV, well below the Gamow peak for this reaction. 
Across the Gamow peak, their measured data can be characterized by 

S{0) = 0.216 ± .006 eV-barns 

S'(0) = 0.0059 ± .0004 barns. 

The rate of this reaction in the solar interior is so much faster than that for 
the pp reaction, that deuterium is effectively instantaneously converted to 
^He, and any uncertainty in this reaction rate has essentially no impact on the 
pp-chain. 

Other d-hurning reactions: The d(d,n)^He and d(d,p)'^H reactions have cross- 
section factors, 5'(0), which are ~10^ times larger than for the d(p,7)'^He 
reaction. However, because the relative deuterium-to-proton abundance ratio 
is ^ 10~^^ in the solar interior, these other reactions do not play a significant 
role. The S-factor for d(^He,p)^He reaction is more than 10^ times larger than 
that for the d(p,7)^He reaction, but this is more than counterbalanced by the 
small •^He/p abundance ratio (< 10^^) and the higher d-|-'^He Coulomb barrier, 
so that this reaction accounts for less than 10~^ of deuterium burning in the 
solar interior {e.g., [33]). 
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4.3 The ^He(^He,2p) ^He reaction 



In a main-sequence star burning hydrogen via the pp-chain, the energy spec- 
trum of the neutrinos emitted from the core is determined by the relative 
rates of the various '^He-burning reactions (and the relative burning rates of 
any subsequent ^Be-burning reactions). This spectrum determines both the 
fraction of the pp-chain energy (26.73 MeV) which is removed from star's en- 
ergy budget (escaping in the form of neutrinos) — and the counting rates for 
various solar neutrino with differing thresholds and response functions. The 
rates of the ^Yic(^Yic,2^fYLe, 3jje(Q;, ^)7Be, 3He(d,p)^He, and 3He(p,e+i/e)^He 
reactions are discussed below. 

Studies conducted in laboratories at the earth's surface of very low energy 
thermonuclear reactions are limited predominantly by background effects of 
cosmic rays in the detectors, leading typically to more than 10 background 
events per hour. Conventional passive or active shielding around the detectors 
can only partially reduce the problem of cosmic-ray background. The best 
solution is to install an accelerator facility in a laboratory deep underground. 
As a pilot project, a 50 kV accelerator facility has been installed in the un- 
derground laboratory at Gran Sasso, where the flux of cosmic-ray muons is 
reduced by a factor 10^ compared with the flux at the surface (Fiorentini 
et al. [34]). This unique project, called LUNA (Laboratory for Underground 
Nuclear Astrophysics) , was designed primarily for a renewed study of the re- 
action ^He(^He,2p)^He at low energies, aiming to reach the solar Gamow peak 
at i?o ± A/2 = 21 ±5 keV. This goal has been reached (Bonetti et al. [35]) 
with a detected reaction rate of about 1 event per month at the lowest en- 
ergy, = 16 keV, with cr = 20 fb or 2 x 10~^^ cm^. Thus, the cross section 
of an important fusion reaction of the hydrogen-burning proton-proton chain 
has been directly measured for the first time at solar thermonuclear energies 
(Fig. 4): extrapolation is no longer needed in this reaction. Across the Gamow 
peak, after correction for electron shielding (see Eq. (6) and its associated 
text) their measured data can be characterized by [35] 

5(0) = 5.32 ± 0.08 MeV-barns 

S'{Q) = -3.7 ± 0.6 barns. 

4.. 4 The ^He (a^^) ^ Be reaction 

The relative rates of the ^lle(^He,2p)^He and ■^lle(a, 7)^Be reactions determine 
the branching ratio between the ppl and ppll+pplll terminations of the pp- 
chain (see Figure 1). The former produces only the low energy p+p— ^d-l-e^-l-i/e 
neutrinos {E,^ < 420 keV), while the latter produce ^Be electron capture neu- 
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trinos and the high-energy neutrinos {E^^^ ~ 14.02 MeV), the neutrinos 
measured in the Homestake, SNO and Super-Kamiokande experiments. 

Measurements of the ^He(Q;, 7)^Be reaction (six via measurements of the direct- 
capture gamma rays and three via measurements of the residual ^Be activity) 
were reviewed by a study group at the INT in 1997 (Section V in Ref. [29]), 
with the result that - 

S{0)DirectCapture^ 0.507 ± .016 kcV-b 

S{0) Residual Activity^ 0.572 ± .026 kcV-b 

SiO)TotalDataSet^ 0.533 ± .013 kcV-b. 

This study group concluded [29] that the grouping of the two data sets (direct 
vs. activity) around two different centroids "suggests the possible presence 
of a systematic error in one or both of the techniques. An approach that 
gives a somewhat more conservative evaluation of the uncertainty is to ... 
determine the weighted mean of those two results (the direct-capture and 
residual-activity data sets). In the absence of information about the source 
and magnitude of the excess systematic error, if any, an arbitrary but standard 
prescription can be adopted in which the uncertainties of the means of the 
two groups (and hence of the overall mean) are increased by a common factor 
of 3.7 (in this case) to make x^=0.46 for one degree of freedom, equivalent 
to making the estimator of the weighted population variance equal to the 
weighted sample variance. The uncertainty in the extrapolation is common to 
all the experiments and is likely to be only a relatively minor contribution to 
the overall uncertainty. The result is our recommended value for an overall 
weighted mean: 

534(0) = 0.53 ± .05 keV-b." 

Now that the ^Bc(p,7)*^B reaction rate seems to be on a firmer footing (see 
Section 4-5 below), the ■^He(a;, 7)'^Be reaction rate, because of the inconsis- 
tency of the direct-capture and ^Be activity results, has become the most 
important nuclear uncertainty in the solar neutrino problem. To try to re- 
solve this issue, at least four new studies of the ^He(Q;, 7)^Be reaction rate are 
currently underway. The approaches include new direct-capture and activity 
measurements, as well as a Coulomb-dissociation experiment. A recent series 
of ^Be-activity measurements [36] at four energies in the range 420 < Ecm < 
950 keV determines an extrapolated intercept of 5*34(0) — 0.53 ± .02 ± .01 
keV-b. 

It should be noted that the ^He(a,7)^Be reaction exhibits a particularly sim- 
ple, non-resonant behavior from its threshold to Ecm ^ 2.5 MeV, and its 
entrance-channel phase shifts are well determined over this energy range. Thus 
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its energy dependence provides a good test of nuclear models. Section V of 
Ref. [29] provides a good review of theoretical work on this reaction. 



4.5 The ^He(p,e^i'e)^He reaction 



The "hep" reaction proceeds through the same weak-interaction as the pp 
reaction (Section above), and as such its cross section factor would be 
expected to be orders of magnitude smaller than for the ^He(^He,2p)^He and 
'^Hc(a,7)^Be reactions and of no interest gnificant termination for the 

pp-chain. However this weak branch produces energetic neutrinos (^E^"-^ = 
18.77 MeV) that extend beyond the endpoint of the neutrino spectrum, 
so that in principle a weak flux could be detected. Both SNO and Super- 
Kamiokande have placed interesting limits on this flux; see Section 5, below. 

The nuclear physics of this reaction is rather subtle: in a naive shell-model de- 
scription, the incident proton would, by the Pauli exclusion principle, occupy 
a state orthogonal to the two s-shell protons in ^He. In the allowed approxima- 
tion, the conversion of the proton to a neutron by the weak interaction would 
then produce a state orthogonal to the nai've s-shell ^He ground state. This 
suggests that smaller components in the nuclear wave function and correc- 
tions such as exchange currents might be unusually important in determining 
the cross section. For this reason theoretical uncertainties in estimating the 
rate are considerably larger than in pp beta-decay, despite modern few-body 
techniques that can be applied to the continuum four-body problem [37] . The 
most recent evaluation of the rate of this weak-interaction reaction [27] gives 
an S-factor of 

Shep{0) = (8.6 ± 1.3) X 10-20 keV b, 

corresponding to a hep termination branch for the pp-chain of 2.5x10""^ and 
a hep neutrino flux of ^hep = 7.9 x 10^ cm'^ in BP04 [9]. 

Other ^He-burning reactions: The S-factor S'(O) for ^Hc(d,p)'^Hc is comparable 
in magnitude to that for '^He(^He,2p)'^He reaction, while the much smaller 
relative abundance of deuterium compared to ^He (d/^He <^ lO^^^) more than 
compensates for the higher ^He-|-^He Coulomb barrier. (Deuterium has a very 
short lifetime ~ 1 sec in the solar core. The ^He lifetime increases sharply with 
decreasing temperature, but is typically ~10^ years in the hotter central core 
of the sun [33].) The net result is that the ^He(d,p)*^He reaction accounts for 
less than 0.1% of ^He-burning in the solar interior. 
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4-6 The ^ Be (e , Li reaction 

The neutral-atom lifetime and the branching ratio (the ratio of the transition 
rate to the first excited state ^Li* to the total ^Li*+^Li^s rate) for the electron- 
capture decay of ^Be is well determined (see Table 7.6 in ref. [38]): 

Ti/2 = 53.22 ± 0.06 days 

uj{^U*)l u;(total) ^ (10.44 ± 0.06) %. 

However, one does not want the neutral-atom lifetime, but the lifetime for 
the conditions of the temperature and density in the solar interior - which 
determine the ionization state and the electron density governing this electron- 
capture decay rate. 

Bahcall and Moeller [39] evaluated the rate for capture of a continuum electron 
as well as the contribution from bound-state capture in Be ions that are not 
totally ionized. While continuum capture dominates, inclusion of bound-state 
capture increases the rate by about 22%. They found 

i?(^Se + e-) = 5.60 X 10-9(p/2)(1 + Xh)T6"'/' 

(18) 

x[l + 0.004(T6-16)]s-^ 

Gruzinov and Bahcall [40] used the density matrix formalism {e.g.^ Feynman 
[41]) in an alternative calculation of the rate. The result from this approach, 
which does not separate bound-state capture from continuum capture, agrees 
with Eq. (18) to better than 1%. 

4.7 The ^Be (p^fB reaction 

While the pp-III branch plays only a very minor role (0.013 %) in solar energy 
production, the decay of the ^B produced in the^Be(p,7) reaction is directly 
responsible for the production of essentially all of the solar neutrinos of in- 
terest to the SNO and Super-Kamiokande detectors. Motivated by interest 
in quantitatively understanding the solar neutrino fluxes measured in these 
detectors, experimentalist have mounted five new direct measurements of the 
^Be(p,7) reaction in the last decade, as well as a number of indirect measure- 
ments in which the rate is deduced either via the strength of a proton-transfer 
reaction or via Coulomb dissociation of a ^B beam. 

Except for the original Kavanagh measurement of the high-energy positrons 
[42] , the direct ^Be(p,7) experiments have measured the reaction cross section 
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by detecting the delayed alpha-particles, emitted following the /3-dccay of 
(ri/2 = 0.77 s) to ^Be. The five modern, low-energy (i?cm < 425 keV) 
measurements [43,44,45,46,47] are consistent with each other, and when each is 
fit with the same theoretical model {e.g., DB94 [48]) they generate a weighted 
average of 21.4 ± 0.5 eV b [47]. However, there is an additional uncertainty 
introduced by the variations between models that can be fit to the data and 
then used to perform the extrapolation to 6*17(0). Junghans et al. [47] compare 
12 model calculations (their Fig. 16) all of which cluster around a central value 
of ~22 ± 0.8 eV b. (See also Jennings et al. [49].) Junghans et al. combine 
these experimental and theoretical uncertainties as 

^17(0) = 21.4 ± 0.5 (exp) ± 0.6 (theory) eV b. 

An important systematic uncertainty in these direct experiments has been 
the thickness and composition of the ^Be targets. To circumvent this issue 
an experiment is currently in its final planning stages at ORNL, which will 
utilize a ^Be beam incident on a hydrogen gas target [50]. (Of course, this will 
introduce other, but different, systematic issues, such as the distribution of 
the hydrogen gas in the windowless target and the charge-state distribution of 
the recoihng ^B nucleus in the recoil separator.) Two other proposed methods 
for avoiding the ^Be-target issues are indirect measurements of the reaction 
rate by (a) proton-transfer reactions or (b) the Coulomb dissociation (CD) of 
accelerated ®B ions: 

(a) Transfer-reaction experiments [51] have been used to extract the asymp- 
totic normalization coefficients (ANCs) from the reactions ^°B(^Bc,*^B)^Be and 
-^^N(^Be,^B)^^C coupled with measurements of the ANCs for the -'^'^C(p,7) and 
^Be(p,7) reactions. The ANCs from these reactions determine the behavior 
of the tail of the ^Be+p wave function in the peripheral region that domi- 
nates capture reactions in such loosely bound nuclei. This method has been 
checked by comparing an ANC measurement of the ^^0(^He,d)^^F reaction 
with direct ^^0(p,7)^^F measurements, yielding agreement at the level of 9% 
[52]. The ANC analysis of the ^0B(7Be,SB)9Be and ^^N(^Bc,«B)i3C reactions 
by Azhari et al. [51] gave ^17(0) = 17.3 ± 1.8 eV b. An ANC analysis [53] of 
a number of ^B break-up studies (involving a variety of energies and targets) 
yielded a value of 5'i7(0) = 17.4±1.5eVb. However, it should be noted that 
a difi^crent, separate analysis [54,55] of the ^^C(^B,^Be)"'^^N break-up reaction 
gave 5'i7(0) = 21.2 ± 1.3 cV b, indicating that substantial uncertainties remain 
in the interpretation of the measurements. 

(b) Several ^B Coulomb-dissociation experiments (in which the break-up 

products are measured at very forward angles, corresponding to very large 
impact parameters and hence very small nuclear contributions) have been 
carried out for ^B beam energies from 52 to 254 MeV/nucleon. The four most 
recent Coulomb dissociation experiments [56,57,58,59] deduce 5'i7(0) values in 
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the range from 17.8 to 20.6 eV b. In the absence of a separate cahbration of 
the CD method, at this time the comparison of the these CD measurements 
with the direct ''Be(p,7) cross section measurements may be more important 
as a check on the level of validity of the CD method for determining (p,7) 
S-factors. Such a check is important in assessing CD for other cases where 
direct measurements may not be possible. 



4.8 The Carbon-Nitrogen Cycle 

Hydrogen burning via the CN-cycle (Fig. 2) accounts for less than 2% of the 
^He produced in the sun. In their analysis of the status of the various nu- 
clear reactions in this cycle, the 1997 INT study group [29] concluded that 
by far the most significant uncertainty in this cycle was the poorly defined 
role of the subthreshold resonance [Ecm = -504 keV) in the extrapolation of 
the measured rate of the ^''N(p,7)^^0 reaction. This reaction is ~100 times 
slower than the other reactions in the CN-cycle and therefore determines the 
overall rate of this cycle. Based on their analysis of the data of Lamb & Hes- 
ter [60] and Schroder et al. [61], the INT workshop participants recommended 
Stot{^) = 3.5^2^0 keV b, and emphasized the need for new experiments to 
improve our understanding of the ^''N(p,7)^^0 reaction. Subsequently, in a 
further R-matrix reanalysis of the Schroder data, Angulo and Descouvement 
[62] extracted a reduced proton width for the subthreshold resonance corre- 
sponding to the 6.79-MeV state that was much smaller than the one used in 
the original Schroder analysis. This smaller reduced proton width resulted in 
a value of Stot{0) = 1.77 ± 0.20 keV b for this reaction [62]. 

Recently, preliminary results have become available from two new, indepen- 
dent measurements of the rates for each of the gamma-ray transitions for the 
i'^N(p,7)i-'50 reaction, at the LENA (UNC-TUNL) facility and at the LUNA 
(Bochum-GranSasso) facility. The LENA measurements [63] utilized an im- 
planted target (^^N atoms implanted in 0.5-mm thick Ta backings) and were 
carried out down to Ecm — 145 keV. The LUNA measurements [64] utihzed 
both a solid target (for Ecm > 130 keV) and a differentially-pumped N2 gas 
target (for Ecm > 70 keV). These LENA and LUNA data sets are in good 
agreement, well within their statistical uncertainties, and when combined with 
a remeasurement [65] of the lifetime of the subthreshold state at E^ = 6.79 
MeV, these new measurements indicate a factor of ~2 reduction in the Stot{0) 
value for this reaction, to 1.67 keV b, compared to the Schroder value of 3.5 
keV b [61,29]. (Preliminary reports for these two data sets indicate statistical 
uncertainties of ~ 5% and systematic uncertainties of ~ 9% in their determi- 
nation of Stoti^)- However, it should be noted that the analysis of the 70-keV 
LUNA data point is not yet complete.) 
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A preliminary examination [63,64] of the consequences of this factor of ~2 
reduction in the rate of the ^^N(p,7)^^0 reaction finds (a) a corresponding 
factor of ~2 reduction in the CN component of the solar neutrino flux and 
(b) an increase of ~1 Gyr in globular cluster ages based on turn-off points 
from the main-sequence. When the data are final, the effects on detailed SSM 
neutrino fiux predictions can be evaluated more quantitatively. 



5 Solar Neutrino Fluxes: Theory and Experiment 

The nuclear cross sections discussed here, together with other solar model 
parameters {e.g., abundances, the solar age), can be incorporated into solar 
model calculations, leading to predictions that can be verified in helioseismol- 
ogy or neutrino fiux measurements. The principal neutrino fiuxes from the pp 
chain are the pp and beta-decay sources and the line sources (broadened 
by about 2 keV from thermal effects) from electron capture on ^Be. As seen 
from Figure 1, one of the reasons for the importance of these fiuxes is that they 
provide a direct test of the competition among the ppl, ppll, and ppIII cycles 
that comprise the pp chain: the ^Be and ^B neutrinos tag the ppll and ppIII 
cycles, respectively, while the pp neutrinos then determine the total rate of 
fusion in the sun. This competition is an important check on the SSM, specif- 
ically its prediction of the sun's central core Tc. For example, the ^B fiux is 
found to vary as ~ T^'^ when solar model input parameters are varied (while 
conserving the luminosity), while the fiux ratio (j)(^Be)/(j){^B) ~ T'^"^. The 
precision with which these fiuxes can be predicted in the SSM is also crucial 
to current efforts to better determine the parameters governing solar neutrino 
oscillations. 

Results from the most recent Bahcall and Pinsonneault SSM (BP04) [9] are 
summarized below, where ElO = lO^*^: 



Source 


BP04 (cm-'^s-'^) 


BP04+ fcm-^s-i) 


PP 


5.94(1±0.01) ElO 


5.99 ElO 


pep 


1.40(1±0.02) E8 


1.42 E8 


hep 


7.88(1±0.16) E3 


8.04 E3 


^Be 


4.86(1±0.12) E9 


4.65 E9 


8B 


5.79(1±0.23) E6 


5.26 E6 


13n 


5.71(1±0.36) E8 


4.06 E8 


150 


5.03(1±0.41) E8 


3.54 E8 


17p 


5.91(1±0.44) E6 


3.97 E6 
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These results reflect the dominance of the ppl cycle in the pp-chain for 
present-day solar conditions: ~ 85% of the sun's energy in generated via 
this path. The ppll (~15%) and ppIII (~0.02%) account for the remainder. 
Approximately 1.7% of ^He synthesis occurs via the CNO cycle. Cross section 
uncertainties lead to substantial ranges in the predicted fluxes. 

The cross section uncertainties discussed in this paper are included in the 
overall model uncertainties quoted for BP04, which uses older element abun- 
dances from Grevesse and Sauval [66]. Also shown are the results for BP04-I-, 
a calculation identical to BP04, but using more recent analyses of the solar 
surface abundances (C, N, O, Ne, Ar) of AUende Prieto et al. [67]. Contribut- 
ing to the uncertainties in extracting solar surface abundances are systematics 
due to blending of atomic lines and uncertainties in modelling the solar atmo- 
sphere. The AUende Prieto et al. abundances substantially lower the surface 
heavy element to hydrogen ratio (from 0.0229 to 0.0176), which leads to a 
shallower convective zone and discrepancies with helioseismology. The BP04-I- 
neutrino fluxes reflect a shift in the pp chain toward the ppl cycle because 
the lower metalicity produces a somewhat cooler sun. Because of its tempera- 
ture sensitivity, the ^B flux shows the largest change among pp fluxes, a 10% 
decrease. The lower abundances for C, N, and O directly influence CN-cycle 
neutrino fluxes, as does the lower value for Tc, leading to changes on the order 
of 40%. BP04 neutrino spectrum is shown in Fig. 5. The ^B neutrino spectrum 
deviates slightly from an allowed shape because the flnal ®Be state is a broad 
resonance. 

The prospect of directly probing the thermonuclear reactions occurring in the 
solar core helped to motivate Ray Davis and his colleagues [68] to mount 
their pioneering chlorine solar neutrino experiment. Sited nearly a mile un- 
derground in the Homcstake Mine to avoid cosmic ray-induced backgrounds, 
this radiochemical detector operated almost continuously from its construc- 
tion in 1967 until 2002. The experiment exploited a reaction, ^^Cl(z/e)C^)'^''Ar, 
that had been first suggested by Pontecorvo [69] and by Alvarez [70]. The 
technique involved the collection of the product '^''Ar (ri/2 ~ 35 d) from a 
tank containing 0.61 kiloton of perchloroethylene (C2CI4). Gas was removed 
from the tank about every two months by a helium purge, then circulated 
through a condenser, a molecular sieve, and a charcoal trap cooled to the 
temperature of liquid nitrogen. The resulting efficiency for collecting the ^^Ar 
was high, typically greater than 95%, as was demonstrated by the recovery 
of known amounts of the carrier gases '^^Ar or '^^Ar that were introduced at 
the start of each run. After extraction, the trap was heated and swept by 
helium. The extracted gas was passed through a hot titanium filter to remove 
reactive gases, and the argon was then separated from other noble gases by 
gas chromatography. The purified argon was loaded into a miniaturized gas 
proportional counter, where the subsequent electron capture on ^^Ar was mea- 
sured through the 2.82 keV of energy that is released as the atomic electrons 
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in ^ CI adjust to fill the K-shell vacancy. The counting typically continued for 
about a year (~ 10 half lives). 

After accounting for small contributions from neutron- and comic-ray-induced 
backgrounds, Davis determined a solar neutrino production rate of ~ 0.5 ^^Ar 
atoms/day in the tank. Despite the low counting rate, the final results for the 
experiment achieved an impressive accuracy, 2.56 ± 0.16 (stat) ± 0.16 (syst) 
SNU, where a Solar Neutrino Unit (SNU) is lO-^^captures/^'^Cl atom/s [71]. 
The chlorine experiment is primarily sensitive to neutrinos (~ 76%), which 
can induce the strong super-allowed transition to the 4.99 MeV state in ^'''Ar, 
and ^Be neutrinos (~16%), which arc sTifficiently energetic to overcome the 
814 keV threshold of the ground-state transition. The discrepancy between 
the chlorine results and the SSM predictions (8.5 ± 1.8 SNU for BP04) was 
the genesis of the solar neutrino problem. 

Two similar radiochemical experiments, SAGE and GALLEX, began solar 
neutrino measurements in January 1990 and May 1991, respectively. These 
exploited the reaction ^^Ga(i/e,e~)^^Ge which, because of the strength and 
low Q-value for the ground-state transition, is primarily sensitive to the low- 
energy pp neutrinos. SAGE, which employed 60 tons of liquid Ga metal, con- 
tinues to operate in the Baksan Neutrino Observatory, under 4700 mwe of 
shielding provided by Mount Andyrchi in the Caucasus. GALLEX and its 
successor GNO, which used 30.3 tons of Ga as GaCls in a hydrochloric acid 
solution, operated until very recently in the Gran Sasso Laboratory in Italy 
(3300 mwe). 

The primary challenge in both experiments is the more complicated chem- 
istry of Gc extraction, which is done about every three weeks. In SAGE the 
Ge is separated by vigorously mixing into the gallium a mixture of hydro- 
gen peroxide and dilute hydrochloric acid. This produces an emulsion, with 
the Ge migrating to the surface of the emulsion droplets where it is oxidized 
and dissolved by hydrochloric acid. The Ge is extracted as GeCl4, purified 
and concentrated, synthesized into GeH4, and further purified by gas chro- 
matography. The overall efficiency, determined by introducing a Ge carrier, is 
typically 80%. In GALLEX/GNO the Ge is recovered as GeCh by bubbhng 
nitrogen through the solution and then scrubbing the gas through a water 
absorber. The Ge is further concentrated and purified, and finally converted 
into GeH4. The overall extraction efficiency is typically 99%. 

In both experiments the GeH4 is inserted into miniaturized gas proportional 
counters, carefully designed for radiopurity, and the ^^Ge is counted as it de- 
cays back to ^^Ga (ri/2 = 11.43 d). As in the case of "^''Ar, the only signal is 
the energy deposited by Auger electrons and X rays that accompany atomic 
rearrangement in Ga. Both K and L captures can be detected. Intense ^^Cr 
neutrino sources were used to calibrate both detectors [72] . SAGE is currently 
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conducting a test calibration with a ^^Ar source, in preparation for an antici- 
pated repeat of this experiment with ~ 2 MCi of ^^Ar [73] . 

The most recent results of SAGE give a counting rate of 66.9 ± 3.9 (stat) ± 
3.6 (syst) SNU [74], while the GALLEX/GNO result is 69.3 ± 4.1 (stat) ± 
3.6 (syst) SNU [75]. The rates arc in good agreement, but well below the SSM 
prediction (131 ± 11 SNU in BP04). 

Two types of direct-counting experiments have also been done, using the 
water Cerenkov detectors Kamiokande and Super-Kamiokande and the heavy- 
water detector SNO (Sudbury Neutrino Observatory). Kamiokande, a 4.5- 
kiloton cylindrical imaging water Cerenkov detector, was originally designed 
for proton decay but was later reinstrumented to detect solar neutrinos. 
(The improvements included sealing the detector against radon inleakage and 
recirculating the water through ion exchange columns.) It operated for nearly 
a decade (1987-95) detecting solar neutrinos by the Cerenkov light produced 
by recoiling electrons in the reaction z/^^ + e — > z/^ + e'. Both electron and 
heavy-flavor neutrinos contribute to the scattering, with aivg) / aivn) ~ 7. The 
threshold of the detector was initially 9.5 MeV but was ultimately lowered to 7 
MeV, due to Kamiokande III improvements in electronics and the addition of 
wavelength shifters to improve PMT light collection. The inner 2.14 kilotons of 
water was viewed by 948 50-cm photomultiplier tubes (PMTs), which provided 
20% photocathode coverage. An additional 123 PMTs viewed the surrounding 
1.5m of water, which served as an anticounter. The fiducial volume employed 
for solar neutrino measurements was the central 0.68 kilotons, the region most 
isolated from the high-energy gamma rays generated in the surrounding rock 
walls of the Kamioka mine. Kamiokande found a *^B neutrino flux of (2.80 ± 
0.19 ± 0.33) X lOVcm^s [76], 48% of the BP04 SSM prediction of (5.79 ± 
1.33) X lOVcm^s. 

The experiment was remarkable in several respects. It was the first to measure 
solar neutrinos in real time. It exploited the sharp forward peaking of the 
scattered electrons, in the direction of the incident neutrino, to extract solar 
neutrino events from a significant but isotropic background. The "pointing" 
of events back to the sun provided the first direct evidence that neutrinos 
originated from the sun. 

Super-Kamiokande, the massive 50-kton successor to Kamiokande, collected 
1496 days of solar neutrino data over its first run, from May 1996 through 
July 2001. With its high PMT coverage (40%), lower threshold (~ 5 MeV), 
and much larger fiducial volume (22.5 ktons), Super-Kamiokande collected 
data at a rate ~ 100 times faster than that of its predecessor, Kamiokande. 
The result was an extraordinary set of precise measurements of the fiux, recoil 
electron spectrum, and seasonal (earth's orbital eccentricity) and day/night 
(earth regeneration of neutrinos) time variations of the solar neutrinos. The 
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experiment placed a significant constraint on the hep flux (< 7.3 X 10^/cm^s, 
a hmit about 10 times the BP04 prediction) and on various exotic neutrino 
properties, such as magnetic moments. Super-Kamiokande's '^B neutrino flux 
result is (2.35 ± 0.19 ± 0.33) x lO^cm^s [77], 41% of the BP04 prediction. 

The SNO detector was built deep in the Creighton ^^9 nickel mine in Sudbury, 
Ontario, for the purpose of determining the flavor content of the ®B solar 
neutrino flux. A central acrylic vessel containing one kiloton of heavy water 
is surrounded by 5m (7 ktons) of light water to shield the inner volume from 
neutrons and gammas. The detector is viewed by 9500 20-cm PMTs, providing 
56% coverage. 

The heavy water allows the experimenters to exploit three reactions with 
varying flavor sensitivities 

t'e + d^P + P + e" (CC: charged current) 

Ux + d ^ u'^ + n + p (NC: neutral current) 

J^x + ^~ ^ ^'x + e~' (ES: elastic scattering) 

As the Gamow-Teller threshold for the CC reaction is concentrated close to 
the p + p threshold of 1 .44 MeV, most of the neutrino energy is transferred 
to the electron. Thus the CC electron spectrum provides a better probe of i/g 
spectral distortions than does the ES electron spectrum. 

The NC reaction, which is observed through the produced neutron, provides 
no spectral information, but does measure the total solar neutrino flux, inde- 
pendent of flavor. The SNO experiment has used two techniques for measuring 
the neutrons. In the initial SNO pure D2O phase the neutron signal was cap- 
ture on deuterium, which produces 6.25 MeV gammas. In a second phase 2.7 
tons of salt was added to the heavy water so that CI would be present to 
enhance the capture, producing 8.6 MeV gammas. The NC and CC events 
can be separated reasonably well because of the modest backward peaking 
(~1 - cos(6')/3) in the angular distribution of the latter. This allowed the ex- 
perimenters to determine the total and electron neutrino fraction of the solar 
neutrino flux. A third phase has begun in which direct neutron detection is 
provided in pure D2O by an array of ^He-flUed proportional counters. 

The ES, which we noted previously measures electron and heavy-flavor neutri- 
nos, but with reduced sensitivity to the latter, is the same reaction exploited by 
Super-Kamiokande. Because the SNO volume is smaller, the SNO ES statistics 
are reduced. However, the lower SNO backgrounds lead to a smaller system- 
atic error, so that the overall the SNO and Super-Kamiokande sensitivities 
are comparable. Because the ES and CC scattered electrons are measured in 
the same detector, the ES reaction provides an important cross check on the 
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consistency of the results from the CC and NC channel. 
The SNO results are presented in Fig. 6. The deduced fluxes are [78] 
= (1.59 ± 0.08 ± 0.07) X lOVcm^s 

(l)Es = (2.21 ± 0.28 ± 0.10) X lOVcm^s 

(l)NC = (5.21 ± 0.27 ± 0.38) x lOVcm^s 

The total flux (NC) is in agreement with the SSM prediction, but the CC 
results indicate that almost 70% of the solar neutrinos arrive on earth as 
heavy-flavor neutrinos. The ES results are in excellent accord with this con- 
clusion, as Fig. 6 shows, as well as with the Super-Kamiokande results. The 
data rule out the hypothesis of no flavor change at greater than 7a. 



6 Conclusions 

This paper has reviewed two problems. The flrst was the 70-year effort to 
understand solar energy generation and stellar evolution as a consequence of 
the thermonuclear reactions occurring in the interiors of stars. It led to col- 
laborations between solar modelers, who constructed codes to describe the 
evolution of our sun and to predict its helioscismology and neutrino produc- 
tion, and laboratory experimentalists, who provided the nuclear reaction data 
essential to the SSM. The second was the nearly 40-year program to measure 
the flux, flavor, and spectrum of solar neutrinos. This effort began with the 
pioneering experiment of Ray Davis, Jr. and his collaborators at Homestake, 
and led to the current generation of sophisticated active detectors, SNO and 
Super-Kamiokande. 

Both programs required great patience — years of effort to reflne the SSM and 
to develop better solar neutrino detection techniques — before it became clear 
that the resolution of the solar neutrino problem had to be new physics. That 
resolution included a new phenomenon, matter enhancement of neutrino os- 
cillations, which altered our view of neutrino oscillations and opened up new 
experimental tests, such as spectral distortions and day-night effects. This 
further motivated the development of robust detectors like SNO and Super- 
Kamiokande. 

The conclusion that was reached after decades of work is quite profound: the 
solar neutrino discrepancy first identified by Davis is due to new particle 
physics, neutrino oscillations. This phenomenon lies outside the minimal stan- 
dard model, requiring massive neutrinos and a nontrivial relationship between 
mass and weak-interaction eigenstates. Analyses that take into account not 
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only the solar neutrino data described above, but also the KamLAND reactor 
neutrino data [79], conclude [80] 

5mii = m^-m? ~ (8.2 ± 0.3 ± 0.9) x 10-^ eV^ 

tan2^i2 - 0.39 ± 0.05 ± 0.15, 

where 9i2 is the neutrino mixing angle and mi and m2 the mass eigenvalues. 
The KamLAND data are very important in further constraining the mass 
difference and thus the vacuum oscillation length. 

This explanation of the solar neutrino puzzle raises new questions. What 
accounts for the scale of neutrino masses, a scale so different from other 
standard-model fermion masses? Why are neutrino mixing angles large? (The 
corresponding atmospheric-neutrino mixing angle 623 is maximal, to within 
measurement uncertainties.) There are many parameters in the neutrino mix- 
ing matrix not yet determined. These include the absolute scale of neutrino 
mass, the hierarchy (where the nearly degenerate pair of neutrinos that par- 
ticipate in solar neutrino oscillations are lighter than or heavier than the third 
neutrino), the size of the third mixing angle ^^13, the charge conjugation prop- 
erties of neutrinos (that is, whether the masses are generated from Majorana 
or Dirac mass terms) , and the size of the CP-violating phases that appear in 
the mixing matrix. The answers to these questions could help point the way 
to important extensions of the standard model. 

One of the most remarkable outcomes of the solar (and atmospheric) neutrino 
problems is that the neutrino mass differences that emerged are in a range ac- 
cessible to terrestrial experiments with accelerator and reactor neutrinos, such 
as KamLAND. Despite the importance of future terrestrial experiments — such 
as double beta decay searches for Majorana masses and long-baseline neutrino 
experiments to look for CP violation — solar neutrinos still have an important 
role to play. For the foreseeable future {e.g., until beta beams are developed), 
solar neutrinos will remain our only intense source of electron neutrinos. Some 
future goals that might be reached by exploiting these neutrinos include [81]: 

* A low-energy solar neutrino measurement to test the consequences of a large 
9 12- One feature of the accepted LMA (Large Mixing Angle) neutrino oscilla- 
tion explanation of current data is that the survival probability for low energy 
neutrinos is substantially higher than that for the electron neutrinos mea- 
sured in SNO's CC channel. Verifying this prediction is important. 

* Future high-statistics pp solar neutrino experiments to improve our knowl- 
edge of 612- Solar neutrino experiments are the primary source of information 
on 5mi2 and 612. While KamLAND has helped to narrow the LMA solution 
region in 6ml2 , neither it nor Borexino will appreciably improve our knowledge 
of the mixing angle. A pp solar neutrino measurement with an uncertainty of 
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3% would constrain 612 with an accuracy comparable to that possible with 
the entire existing set of solar neutrino data. Thus a goal of next-generation 
pp solar neutrino experiments is to achieve 1% uncertainty, thus substantially 
tightening the constraints on 9i2- 

* Precise flux measurements of the low-energy (pp an<f Be) solar neutrinos to 
improve limits on sterile neutrinos coupling to the electron neutrino. Sterile 
neutrinos with even small couplings to active species can have profound cos- 
mological effects. The coupling of to sterile states can be hmited by CC 
and NC measurements with an accurately known neutrino source. KamLAND 
(together with the SNO CC and NC data) should ultimately limit the sterile 
component of neutrinos to ~ 13%. This bound can be improved by mea- 
suring the NC and CC interactions of pp and ''Be neutrinos to accuracies of 
a few percent. In general one expects the sterile component of solar neutrino 
fluxes to be energy dependent. Thus low-energy solar neutrino experiments 
are an important part of such searches. 

* Preliminary estimates indicate that limits on the neutrino magnetic moment 
could be improved by an order of magnitude in future pp neutrino ES exper- 
iments. A neutrino magnetic moment will generate an electromagnetic con- 
tribution to neutrino ES, distorting the spectrum of scattered electrons. The 
effect, relative to the usual weak amplitude, is larger at lower energies, making 
a high-precision pp solar neutrino experiment an attractive testing ground for 
the magnetic moment. The strongest existing laboratory limit, //.(t'e) < 1-5 x 
lO"^*^ Bohr magnetons, could be improved by an order of magnitude in a 1% 
experiment. 

* Future pp solar neutrino experiments to probe CPT violation with an order of 
magnitude more sensitivity than has been achieved in the neutral kaon system. 
As neutrinos are chargeless, they provide an important testing ground for 
CPT violation. If CPT is violated, the neutrino mass scale and mass splittings 
will differ between neutrinos and antineutrinos. (This has been offered as a 
possible explanation of the LSND neutrino oscillation results.) A high precision 
measurement of the pp i/gSurvival probability is an important component of 
CPT violation searches. In some models this will test CPT violation at a scale 
< 10~^° GeV, which can be compared to the present bound from the neutral 
kaon system, < 4.4 x 10"^^ GeV. 

* Solar neutrino detectors are superb supernova detectors. Many of the most 
interesting features in the supernova neutrino "light curve" are flavor spe- 
cific and occur at late times, 10 or more seconds after core collapse. Because 
of their low thresholds, flavor specificity, large masses, and low backgrounds, 
solar neutrino detectors are ideal for following the neutrino emission out to 
late times. Solar neutrino detectors will likely be the only detectors capable 
of isolating the supernova i/g fiux during the next supernova. The 3 ms delep- 
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tonization burst, important in kinematic tests of neutrino mass, is mostly of 
this flavor. Electron-flavor neutrinos also control the isospin of the nucleon 
gas — the so called "hot bubble" — that expands off the neutron star. The hot 
bubble is the hkely site of the r-process. Currently we lack a robust model 
of supernova explosions. The open questions include the nature of the ex- 
plosion mechanism, the possibility of kaon or other phase transitions in the 
high density protoncutron star matter, the effects of mixed phases on neutrino 
opacities and cooling, possible signatures of such phenomena in the neutrino 
light curve, and signals for black hole formation. Supernovae are ideal labora- 
tories for neutrino oscillation studies. One expects an MSW crossing governed 
by ^13, opening opportunities to probe this unknown mixing angle. The MSW 
potential of a supernova is different from any we have explored thus far be- 
cause of neutrino-neutrino scattering contributions. Oscillation effects can be 
unravelled because the different neutrino flavors have somewhat different av- 
erage temperatures. All of this makes studies of supernova neutrino arrival 
times, energy and time spectra, and flavor composition critically important. 
Finally, the detection of the neutrino burst from a galactic supernova will pro- 
vide an early warning to optical astronomers: the shock wave takes from hours 
to a day to reach the star's surface. 

* Future solar neutrino experiments could provide crucial tests of the SSM. 
The delicate competition between the ppl, ppll, and pplll cycles comprising 
the pp chain is sensitive to many details of solar physics, including the core 
temperature, the sun's radial temperature profile, the opacity, and the metal- 
icity. We have noted that this competition can be probed experimentally by 
measuring the ''Be flux (which tags the ppll cycle), the flux (which tags 
the pplll cycle), and the pp flux (which effectively tags the sum of ppl, ppll, 
and pplll) . Precise measurements of the pp and ''Be neutrino fluxes, combined 
with further improvements in nuclear cross section determinations, could sig- 
nificantly improve our understanding of current conditions in the solar core. 

* The CN neutrino flux is an important test of stellar evolution. The CN cy- 
cle is thought to control the early evolution of the sun, as out-of-equilibrium 
burning of C, N, and O powers an initial convective solar stage, thought to 
last about 10^ years. Furthermore, one of the key SSM assumptions equates 
the initial core metalicity to today's surface abundances. A measurement of 
CN cycle neutrinos would quantitatively test this assumption. The recent con- 
troversy about surface elemental abundances provides additional motivation. 

This work was supported in part by the US Department of Energy under grants 
DE-FG02-00ER-41132, DE-FC02-01ER-41187 (SciDAC), and DE-FG02-91ER- 
40609. 
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proton-proton chain 



?H+ e- + !H «D + Vj 1.442 MeV E* = 1.442 MiV 

'H+'K-»«D + e+ + v, 1.442 MeV B. 0.429 IteV 

'H -* Me+y 5.493 MeV 

(a) 3He + 3He^^He+ap: l&»M«ii¥ 
:(t>)^He + 'He -+ ■^e + Y 1 .587 Me V 

pp-n (bl) ^e + B-^'Li + v, 0.862 MeV E,= 

■^i+;l^Hi*2*He 17;-347MeF 



Fig. 1. The pp-chain, showing the three terminations: ppl, ppll, and ppIII. 




Fig. 2. The primary CN-cycle, together with three subsidiary cycles which can 
bypass the ^^N(p,a)^^C termination to extend this processing as far as ^'^Ne. 
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Fig. 3. The energy dependence of the integrand in Eq. (6), defining the parameters 
of the Gamow peak. 
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Fig. 4. The measured energy dependence of the astrophysical S(E) factor for the 
^He(^He,2p)*^He reaction [35]. The sohd curve is a fit to the data including the 
electron shielding parameters deduced by Aliotta et al. [20] from a study of the 
^He(d,p)^He reaction; the dotted curve represents S{E) for the bare nuclear in- 
teraction, extracted from the solid curve by removing the effects of the electron 
shielding. 
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1 2 3 4 5 6 

Fig. 6. The flux of solar neutrinos which have either or r flavor vs. the corre- 
sponding flux of electron neutrinos, as deduced from the charge-current (CC), neu- 
tral-current (NC), and elastic-scattering (ES) reactions in the SNO heavy-water 
detector [78]. 
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